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Abstract 



We study the structure and stability of discrete breathers (both pinned and mo- 
' bile) in two-dimensional nonlinear anisotropic Schrodinger lattices. Starting from a 

set of identical one-dimensional systems we develop the continuation of the localized 
pulses from the weakly coupled regime (strongly anisotropic) to the homogeneous 
^ ' one (isotropic). Mobile discrete breathers are seen to be a superposition of a localized 

■ mobile core and an extended background of two-dimensional nonlinear plane waves. 

^ , This structure is in agreement with previous results on one-dimensional breather 

^ I mobility. The study of the stability of both pinned and mobile solutions is performed 

using standard Floquet analysis. Regimes of quasi-collapse are found for both types 
of solutions, while another kind of instability (responsible for the discrete breather 
fission) is found for mobile solutions. The development of such instabilities is studied 
examining typical trajectories on the unstable nonlinear manifold. 
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1 Introduction 



We present a numerical study of exact breather dynamics in two-dimensional 
nonlinear Schrodinger lattices. This issue fits well with the enduring scien- 
tific interests of Serge Aubry in the many faceted subject of localization and 
transport in nonlinear macroscopic discrete systems, where Serge's outstand- 
ing contributions are widely recognized. Many ideas and lines of study followed 
in the investigation reported here have found a source of inspiration in early 
works of Aubry and collaborators, among which the doctoral works of Th. 
Cretegny and J.L. Marm deserve special mention. 

The existence and properties of localized solutions in extended discrete sys- 
tems have atracted interest in a broad range of physical fields [1]. Discrete 
breathers, sometimes referred to as intrinsic localized modes, are spatially lo- 
calized and time periodic sohitions. These solutions arise in the context of 
nonlinear discrete systems and are of fundamental interest for varied physi- 
cal applications such as pulse propagation in nonlinear optics, energy storing 
and transport in biomolecules, plasma physics, etc... The existence of dis- 
crete breathers in these systems has been proved rigorously [2] for a number 
of equations with physical relevance and, contrary to continuous nonlinear 
equations, their existence can be regarded as a generic feature of these sys- 
tems. One of the most important class of equations are the so-called discrete 
nonlinear Schrodinger lattices [3,4]. The existence of discrete breathers has 
been proven for a wide range of systems belonging to this class of nonlinear 
difference-differential equations. In particular, the most important example of 
wide aplicability is the standard nonlinear Schrodinger equation. For instance, 
this equation was employed in [6,7] for describing the propagation of localized 
beams in an array of nonlinear (Kerr type) waveguides, having experimental 
validation subsequently reported in [8,9]. 

The study of two-dimensional nonlinear Schrodinger lattices has atracted 
much atention [10,11] in recent years due to the new phenomena emerg- 
ing when the dimensionality of the lattice is increased. Some examples of 
these new features are the existence of vortex-breathers [12] which supports 
energy flux, the appearance of an energy threshold for the creation of dis- 
crete breathers [13,14,15,16,17] and the ubiquity of an instability (the quasi- 
collapse) of some discrete breather solutions leading to a highly localized pul- 
son state [18,19,20,7,22,23]. These theoretical efforts have their counterpart in 
recent advances in the field of nonlinear optics. The studies of two-dimensional 
arrays of coupled nonlinear waveguides allow the experimental observation of 
those effects studied theoretically. Specially relevant is the recent experimental 
breakthrough (theoretically designed in [24]) by Fleisher et al [25,26], where a 
two-dimensional array of nonlinear waveguides is induced in a photosensitive 
material. This technique provides a clear experimental verification of the two- 
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dimensional discrete breather existence in this system. In particular, besides 
the observation of standard discrete breathers, these works reported the first 
observations of staggered discrete breathers. 

Our study here will focus on the computation of numerically exact dicrete 
breathers in two-dimensional anisotropic nonlinear Schrodinger lattices, i.e. 
where the couplings in the two spatial directions are different. The use of 
shooting methods allow us to find these solutions and analyze their structural 
and stability properties. Both pinned and mobile discrete breathers are stud- 
ied. In the latter case we will study only the ones whose motion is along one 
axis of the lattice. The analysis of the numerically exact solutions shed light 
on some features of the properties and stability of localized solutions reported 
in previous works. 

The plan of the paper is as follows: In section 2 we provide the technical 

background and details needed for self-contained purposes. First we summa- 
rize in 2.1 the main conclusions on the dynamics of ID Schrodinger discrete 
breathers reported in [27,28]. A detailed account of the numerical methods 
(SVD-regularized Newton continuation of operator fixed points) that we have 
used can be found in that reference. Also in 2.1, we discuss briefly the most 
relevant formal differences with respect to alternative approaches to (one- 
dimensional) exact mobility of discrete breathers, e.g. those in refs [29] and 
[30,31]. In subsection 2.2 we introduce the two-dimensional anisotropic Salerno 
lattice and provide explanations on the implementation of the numerical pro- 
cedures used to study the dynamics of 2D discrete breathers. 

The analysis of the results of our numerics on pinned discrete breathers for 
anisotropic nonlinear Schrodinger lattices is reported in section 3. We present 
the numerical computations of the fixed point norm, as a function of three 
parameters: breather frequency, transversal coupling, and nonlinearity (see 
below). They show, as anticipated, the so-called quasi- collapse transition, as- 
sociated with the (well-known) existence of thresholds for the breather norm 
in two-dimensional lattices. We present numerically computed sectors of the 
bifurcation surface. We conclude this section with a brief look at the nonlinear 
dynamics on the unstable manifold, whose typical trajectories have been called 
pulson states. Early numerical work on the 2D quasi-coUapse phenomena in 
isotropic lattices was reported in [22,23] and [20]. A two-year-old account of 
the " state of knowledge" on 2D Schrodinger lattices can be found in Section 
six of [3]. 

In section 4 we show results on a type of mobile breather, namely those mov- 
ing along the direction of stronger lattice coupling constant. The structure of 
each of these mobile exact discrete breathers is that of a localized moving core 
superimpossed on a specific extended state of resonant small amplitude radia- 
tion, the background. We present here the results of an extensive Floquet sta- 
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bility analysis of this type of solutions in two sectors of the three-dimensional 
parameter space, which clearly show the existence of two different transitions. 
The tangent space eigenvectors associated to each of the transitions are pre- 
sented, and the relation of the unstable manifold trajectories to pulson states 
is analyzed afterwards. 

We conclude with section 5, where we briefly review the results obtained and 
illustrate their possible implications for mobility of 2D discrete breathers. 



2 Numerical continuation and preliminary results 



The use of numerical tools for the continuation of discrete breather solutions 
has been widely employed since their existence proof was reported (see e.g 
[32,33]). In particular, the design of numerical techniques for finding exact 
mobile breathers based on those employed for the pinned ones has been ex- 
plored in the recent years [27,28,34,35,36] and paves the way to resolving the 
still open question about discrete breather mobility. Here, after reviewing the 
most important features of ID mobile breathers in nonlinear Schrodinger lat- 
tices, we briefly explain how the continuation method is implemented in our 
system. 



2.1 Mobility of one- dimensional nonlinear Schrodinger discrete breathers 



As was introduced in section 1, exact mobility of dicrete breathers in ID 
Schrodinger nonlinear lattices were numerically studied by the authors in pre- 
vious works [27,28]. In particular, a numerical continuation from the integrable 
Ablowitz-Ladik lattice, A-L, (where exact mobile breathers can be calculated 
analytically [37]), 



i$„ = -(^n+l + ^n-l] 



2' ' . 



(1) 



to the standard discrete nonlinear Schrodinger equation, DNLS, 

i6„ = -($„+! + $„_i)-7|$„p$„, (2) 



was performed within the so-called Salerno model 



i$„ = -($n+i + ^n-i) 1 + A*|^ 



(3) 
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where 7, // and v are the parameters accounting for the strength of the non- 
hnear terms. The above equation includes the former two relevant equations 
(1) and (2), for (/^ = 7/2, i/ = 0) and (/^ = 0, = 7/2) respectively, pro- 
viding the desired interpolation needed to develop the continuation scheme. 
The integrability of equation (1) provides for 7 > a two-parametric family 
of discrete breathers 

= sinh/3 sech[/3(n - x{t))] exp[i(Q;(n - x{t)) - n{t))] , (4) 



where the two parameters a e [— tt : tt] and /3 > describe the velocity and 
internal frequency of the solution 



Vb = x = 2 sinh/? sina/ /? , (5) 
LOb — O, — 2 cosh/3 cosa -|- avb ■ (6) 

The Salerno lattice (3), possess two dynamical invariants, namely the hamil- 
tonian. 



n A* n 

+ 24EM1 + A^l*nr), (7) 

A* n 

where denotes the complex conjugate of and the norm 

A/'=-X:Ml + /^|$n|'). (8) 
A* n 

In order to find exact mobile discrete breathers in the Salerno model we define 
a (p, g')-resonant solution referred to some time scale r such that 

K{to) = $n+p(to + qr) . (9) 



Within the above definition, a mobile breather that translates p sites after q 
periods of the internal oscillation will satisfy equation (9) when r = Tj,. In 
this sense, the continuation focuses on the families of {p, g)-resonant discrete 
breathers, that is breather solutions with the two characteristic time scales 
(corresponding to the breather velocity, v^, and the internal frequency, Ub) 
being commensurate. In all the computations wc have used finite lattices with 
periodic boundary conditions (PBC) so that ^n+i = and $o = ^Af (with 
N being the lattice size). 
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In the previous works [27,28] the authors start from those A-L sohitions which 
are [p, g)-resonant, we discrctizc (fine grid) the path /x + z/ = 1 (with fi and 
being positive) along the Salerno model, and go through it computing the 
corresponding (p, g)-resonant solutions for the pairs (/x, u). We can choose the 
path without loss of generality because of the scaling property of equation (3). 
Then, each solution is numerically computed as a fixed point of the map 

M = L^%,,.^ , (10) 



where L is the lattice translation operator L({$„(to)}) = {^n+i{to)}, and 
T{ij,-,^v) is the Tfj-evolution map (T5 = 27i/uJh) following the dynamics dic- 
tated by equation (3) for the corresponding value of (/x = 1 — u)] i.e., 

T(.„.)[{<f„(to)}] ={*n(t0 + T5)}. 

The continuation was then performed for a fine grid of frequencies belonging to 
the family oi {p = l,q = l)-resonant discrete breathers. The most important 
conclusion about discrete breathers in these nonlinear Schrodinger lattices is 
that mobility in the non-integrable regime {i/ ^ 1) demands the presence of an 
extended background to which the fixed point solution is spatially asymptotic 
{n — > 00), i.e. the solution is exactly written as 

*n(t) = *r"W + *I?'^SW- (11) 



This expression defines the purely localized component <l>^'^^^(t) of the solu- 
tion. The background is a finite linear combination of nonlinear plane waves, 

$n (t) = Ylj=i ^j^xp[i(kjn — uj{kj, Aj)t)]. These plane waves are exact 
solutions of the Salerno model (3) being the nonlinear dispersion relation 
u;{kj, Aj) — — 2[1 -I- A^l^j P] cos kj — 2i/\Aj\'^. The results concerning the charac- 
terization of the background are discussed in detail in [27,28], here we briefly 
summarize the most important features: 

( i ) The set of "s" plane waves which take part in the background of a (p, g) —resonant 
discrete breather with internal frequency cub is derived by the simple selec- 
tion rule for the wave-numbers kj 



(jjb q \27r 

i.e. only the plane waves which are {p, g)-resonant with the internal period 

of the breather can be components of {$n (^)}- The number of solutions 

of (12) fixes "s". 

(a) The amplitudes {Aj} of the nonlinear plane waves differ by orders of mag- 
nitude. 
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(in) There exist a strong positive correlation between the amphtude of the back- 
ground and the strength of the Peierls-Nabarro barrier arising from the pe- 
riodic lattice. This correlation is particularly clear when symmetry breaking 
transitions occur for the also studied case oi u < and /x > (see reference 
[28]), and reflects the hnk between non-integrabihty and the existence of 
the background dressing of the mobile core. 

(iv) Finally, the interpretation of the correlation described in (in) is reinforced 
from a study of the energy evolution of the mobile core: There is an energy 
balance brought by the background when the core moves along the lattice. 
In particular, it can be observed how the core energy oscillates periodically 
so that it takes the maximum energy value when the core visits the inter- 
site configuration. This extra energy periodically obtained by the core is 
provided by the background, with the energy maximum clearly related to 
the background amplitude. 

Before concluding this subsection, it is worth commenting on some of the dif- 
ferences between the Newton continuation of fixed points that we use in this 
paper, and other important recent approaches to breather numerics. The work 
by Ablowitz et al [29] uses discrete Fourier analysis to obtain a nonlinear non- 
local integral equation, from where the " ... soliton is thus viewed as a fixed 
point of a nonlinear functional" (sic) in the Fourier transformed space of func- 
tions. Following these authors, their results seem to differ from those of early 
pioneering work [39] (nowadays textbook material [1]) "in which a continuous 
travelling solitary waves were reported using Fourier scries expansions with 
finite period L while assuming convergence as L — > oo" (sic). Ablowitz et al 
term continuous a solution that can be defined off the lattice points, which 
they see as "necessary when discussing travelling waves in lattices" (sic), and 
disagree with some conclusions reported in the earlier works. 

The (" orthodoxy matters" ) discussion above helps us to clarify how differently 
our numerical approaches "sees" the discrete Schrodingcr breather problem: 
The very concept of a variable defined off the lattice points is intrinsically alien 
to our discrete approach, which neither needs of it nor excludes its eventual 
consideration. In contrast to those views (but not at all in logical opposition), 
we consistently view the thermodynamical limit {N — > oo) in lattice space, 
much in the sense used e.g. by Serge Aubry in his celebrated work on the 
Frenkel-Kontorova ground state problem [40]: The infinite size limit is built 
up from a subsequence of PBC (finite) lattices for which the limit is well 
defined. This will make the Fourier-transformed /c-space continuum. 

Closer to our approach in some respects, though technically different in many 
others, is the formal approach purposed recently by James and collabora- 
tors [30,31]. This approach, which uses recent central manifold theorems, was 
brought to our attention very recently and unfortunately, we have to defer to 
a future publication a comparison of our ID computations with theirs, this 
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is beyond the scope of this paper, namely, the anisotropic 2D Schrodinger 
lattices, which we introduce in the next subsection. 



2.2 Two-dimensional anisotropic lattices. 

Motivated by the method described above, we can extend the continuation 
scheme for calculating exact discrete breathers in higher dimensional systems. 
In particular we focus on the two-dimensional nonlinear Schrodinger lattice 

i^nm = - Ci($„+i,^ + $n-l,m)(l + l^\'^n,m\'^) 

n,m\ ) 

-2i/$„,^|$„,^P . (13) 

This lattice can be viewed as the two-dimensional Salerno model. The two 
coupling parameters Ci and C2 provide a technical advantage for numerics 
(see below), but they are also introduced for theoretical and experimental 
interest. The possibility of controlling the ratio between the two linear cou- 
plings of the two transversal directions has been studied in various works as a 
way of analysing how the intrinsic 2D phenomena (such as the quasi-coUapse) 
emerge. In fact, for Ci << C2, = and v — 7/2 equation (13) describes 
a set of weakly coupled nonlinear wavequide arrays and can be considered 
of ^intermediate dimensionality'\ This extreme has been studied 
experimentally in [41] and using perturbative methods in [42]. On the other 
hand, this equation incorporates, as two particular limits, the physically rel- 
evant standard two-dimensional DNLS equation {/i — 0, u — 7/2) and the 
two-dimensional A-L lattice (/x = 7/2, u = 0). The continuation between 
these two limits provides a useful tool for studying the interplay between the 
on-site and inter-site nonlinearity in the 2D anisotropic lattice. Moreover, the 
anisotropy (or freedom in the values of the coupling parameters Ci and C2) 
allows to include an integrable model among the members of the family of 
nonlinear lattices described by (13). That is, for v — 0, Ci = and Cj 7^ one 
obtains a set of integrable A-L chains. In this way, every 2D model included in 
(13) is connected with an integrable model where analytic discrete breathers 
are available. In what follows we set 7 = 2 and u — 1 — /i with 1/ > and 
> 0. 

The dynamics (13) can be derived from the Salerno Poisson structure 
m dA dB dA a-B \ / , ,^ .^n 
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with the hamiltonian 



n,m 

n,m+l 



n.m 



- — ^|$„,^|2 + ^^ln(l + /x|$„,„|2) . (15) 

A* n,m ^^ n,m 

As in the ID Salerno model there is also a second conserved quantity, the 
norm, due to the phase invariance of the equations of motion (13) 

U^-Y.\^{l+ll\^n,m?) . (16) 



1^ 



n,m 



In the same manner as in the ID case we will focus on a special set of 2D 
discrete breathers. For this, we have to generalize the definition (9) of a res- 
onant solution in the two-dimensional case. In this context discrete breathers 
solutions are characterized by three time scales. Namely, one associated with 
the internal oscillation ujh and the other two derived from the translation of 
the localization center, i.e. its velocity Vh = (v^jVy). The subset of 3-tuples 
{uih,Vh) that fulfil the (pa;,^^^, g) -resonance condition 



27r Pc, 

Vx — = — (17) 

Ub q 

27r p„ , , 

Vy-^"^, (18) 

(where p^, Py and q are integers) denote the breather solution that can be 
obtained with our continuation method. These solutions are those that after 
q periods of the internal frequency, ^{to + qTb), translates Px and Py lattice 
sites in the x and y direction of the square lattice, respectively, i.e. 

^n,m{to) = ^n+p,,m+py{to + qTj,) , (19) 



where again PBC are applied $7v,+i,m = *i,m, *o,m = '^'iv,,m, ^n,Ny+i = $n,i 
and ^nfi — ^n,Ny (with and Ny being the lattice size in the x and y direction 
respectively). Consequently, a ipx,Py, g) -resonant state will be a solution of the 
following set of equations 

i^(p.,p,,,,c.,,.,Ci,C.)[{^n,m(io)}] =i^^^i^S^r(^,,,,,Ci,C.)[{^n,rn(io)}] = 

= {<^'n,m(t)}, (20) 
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where the operators Lj are the lattice translation in the i-direction, 

La^[{^n,m(to)}] = {^n+l,m(to)} , (21) 
Ly[{^n,m{to)}] = {^n,m+l(io)} ■ (22) 

Besides, T(^uj,,,u,Ci,C2) is the time evolution operator given by equation (13) over 
one period Ti, — 2Tr / cub, 

Ti.„.,CuC,)[{<^n,mm] = {<^n,m{t0 + n)} ■ (23) 

In order to illustrate the 2D time scales resonance we consider the plane wave 
solutions of equation (13): $n,m(^) = Aexp[i{kxn + kym — u!t)]. These solutions 
possess the following nonlinear dispersion relation 

u;{k, A) = 2(Ci cos A;^ + C2 cosky){l + fxA^) - 2uA^ . (24) 

Hence, we can obtain the subset of plane waves which are {px,Py, g)-resonant 
with some time scale r [i.e. after a time qr they have translated and py 
sites in the x and y direction, respectively). Each member of these subsets will 
be labelled by the pair k = {kx, ky) and from the condition (20) it follows that 

— * 

the corresponding set of values of k for each family will satisfy the relation 
.ft^) = i(p.«-^), (25) 

where m is an integer and p ={px, Py) - In figure 1 the corresponding values of 

— * 

k are represented for two resonances of the type {px = l,Py = 0,q = 1) and 
(Px = 'i-,Py = i,g = 1)- 

The method used for solving equation (20) for each resonant 3-tuple {uJh,Vh) is 
the same as in the ID case, which is extensively described in [28]. The implicit 
function theorem assures that a fixed point solution of a map (20) given by 

— * 

C —{Px, Py, Q, <^6, ^, Ci, C2) can be obtained provided that (i) the Jacobian of 
the operator F^[{$„_^(to)}] — I is invertible and (ii)we know a fixed point of a 

map corresponding to an infinitesimally close set of parameters, ^ — 5^ = {px, 
Py, q, ujh — SuJi), V — 5v, Ci — 5Ci, C2 — 5C2). The first demand can be satisfied 
using a singular value decomposition (SVD) [43,44,45] of the Jacobian in order 
to obtain the pseudo-inverse operator. On the other hand, when the second 
condition is fulfilled convergence of the Newton-Raphson iterative scheme is 
guaranteed. For this, we start with a sufficiently good trial solution, {^°,^(io)} 
and solve the equation 

{^KM} = -DF^[{K^M}\-' ■ Fi[{<^^^^)}\ ' (26) 
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-3-2-10123-3-2-10123 

kx kx 

Fig. 1. Wave numbers, k = (kx,ky), of the (1,0,1) (a) and (1,1,1) (b) resonant 
plane waves for m = (see equation (25)). Different values of C2, while Ci is fixed 
(Ci = 1), are shown. The reference time scale for the resonance is set to r = 2.4315 
(w = 2.584). 

in order to obtain {$i^,m(to)} = {$n,m(^o)} + {^^n,m{to)}- We iterate this 
calculations to the desired convergence, and then the solution {$n,m(^o)}; is 
obtained. In our numerics this is the case when 

Fi[{K,M}] < N . 10-'' , (27) 



(where N is the total number of sites in the square lattice) is fulfilled. Once 
the solution is found we use it as the following trial solution, ^(to)}, for 
solving the map (20) corresponding to the next set of parameters ^ ' = ^ + (5^. 

As an additional benefit, this method provides the linear stability of the com- 
puted solution. For this, we only have to compute the eigenvalues of the 
Jacobian of the operator around the solution DF^{^n,m{to)}]- In fact, this 
Jacobian is the extended Floquet matrix of the solution. Then, a solution is 
linearly stable if all the eigenvalues of the corresponding Floquet matrix are 
inside the unit circle. Moreover, the symplectic character of the dynamics (13) 
implies that all the eigenvalues appear in quadruplets (A, A, 1/A, 1/A) and 
thus for a hnearly stable solutions all the eigenvalues of its extended Floquet 
matrix lie on the unit circle. 

There are two possible paths for developing the continuation method depend- 
ing on the choice of the starting point of the continuation. One possibility 
is to start from the full anti-continuum limit, Ci — C2 — 0, where a pinned 
breather solution of frequency cUb is written as 

^n,m{t) = Sn,noSm,mo\ ^ ex.-p{iu}bt) . (28) 
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Starting from the above solution, we can perform the continuation increasing 
the parameters Ci and C2 as usual, and so obtain the whole family of {px = 
0,Py = 0,q = l)-resonant discrete breathers. An alternative path starts from 
the one-dimensional limit, C2 — 0. The choice of this second limit (which 
implies taking as the very initial trial solution of the continuation the whole 
set of ID solutions described in section 2.1) is justified when seeking mobile 
solutions. As stated above, this limit offers the possibility of studying strongly 
anisotropic lattices as a controlled interpolating situation between one and two 
dimensions. On the other hand, employing this strategy we can only obtain 
those solutions which are {p^ — P,Py — 0, g)-resonant, i.e. the two-dimensional 
continuation of those one- dimensional {p = Px, 5)-resonant discrete breathers. 
Hence, the solution from which we start is 

$n,m(t) = (t) , (29) 

where ^n^{t) is the corresponding (p = p^., g)-resonant one- dimensional solu- 
tion. 

In what follows we will employ both continuation paths when we study the 
case of pinned breathers (section 3), and we will show that the results obtained 
are the same when approaching the same limit (the standard two-dimensional 
DNLS). 



3 Two-dimensional pinned discrete breathers 

As we have discussed, we can choose two different starting points for the 
continuation of (0, 0, 1) -resonant fixed points (pinned breathers) of equation 
(20): (i) the full anti-continuum (AC) limit (Ci — C2 — 0), or (ii) the (one- 
dimensional, ID) limit of uncoupled chains (Ci 7^ 0,6*2 = 0), where they can 
be obtained from continuation along increasing values of the parameter u from 
the ID A-L lattice (1). As a test for our codes, we have checked that both 
paths arrive to the same solution. In fact, unique continuations can proceed 
along any path on the plane of parameters (C2, ly) that we have explored. 

Early works [18,19,20] on the isotropic two-dimensional standard DNLS equa- 
tion analyzed the so-called quasi-collapse instahilitj of pinned discrete breathers, 
i.e. the condensation of all the energy into a few modes in discrete nonlinear 
systems, which corresponds to the onset of a singularity (wave collapse) [21] 
in multidimensional continuum models. Subsequent numerical works [22] ex- 
tended these studies to the isotropic 2D Salerno lattice and addressed the 
question of how the instability is affected by the presence of impurity lattice 
sites. 
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As expected, our results further corroborate the existence of quasi-coUapse 
instabihtics in the anisotropic case: The phase diagram in parameter space 
{cUb, C2, I') consists of two regions (stable and unstable) separated by the 
surface of transition. As we perform the continuation of breather solutions 
across the parameter space we scan the Floquet stabihty of the computed 
solution. In figure 2 we present the two stability transition curves in the plane 
{ujh-, C2, v = 1), i.e. the function C^^{uJi,)^ corresponding to the two different 
continuation starts. The continuation from the AC limit is made through 
the path Ci = C2 and the one from the ID limit is made at Ci = 1. The 
convergence of the two paths at C2 = 1 is clearly seen. 

1 

0.8 
0.6 

^2 0.4 
0.2 


1 2 3 4 5 

COb 



C,=1 



Fig. 2. Evolution of the threshold value of the coupling parameter, C^^, as a function 

of the frequency, Wf,, for two different continuations starts. The values of limit the 
region where pinned discrete breathers are linearly stable (unstable for C2 > C^). 
The instability yields a hyper-localized state (quasi-collapse). The continuation from 
the fully uncoupled limit (Ci = C2 = 0) (filled circles) is performed using the path 
Ci = C2. For the continuation (bold circles) from the 1-dimensional limit {C\ = 1, 
C2 = 0) the coupling in the new direction C2 is progressively increased. 



The criterion for stability of the pinned discrete breather solution derived in 
[19,20], 

'^•^^ >0. (30) 



is of a very general character and our numerics illustrate it clearly. On the 
other hand, the Floquet stability analysis detects the dimensionality (and a 
basis in tangent space) of the unstable linear manifold associated with the 
quasi-collapse instability that these exact discrete breathers experience for 
some parameter values. We have computed numerically, for a fine grid of a;;, 
values and a coarser grid of Ci and the function MiuJb, Ciji'), from which 
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we show some sectors in figures 3 and 4. 

In figure 3 we show the numerically computed norm (16) as a function of the 
breather frequency Af{ujh), for three different values of the transversal cou- 
pling C2, and a fixed value of u = 1 (anisotropic DNLS limit). We observe 
the existence of a minimum value, min J\f{uJb) = 7^ 0, which is thus seen 
as an excitation threshold for the creation of these solutions. The position of 
the minimum ujf^{C2), which naturally increases with C2, separates the stable 
and unstable branches of pinned breathers. Breathers corresponding to values 
of oji) where Af{uJb) has a negative slope are unstable: This is shown in the 
insets, where the Floquet spectra of two representative examples of pinned 
discrete Schrodinger breathers are plotted in the complex plane. Note that 
the high accuracy of the numerical solution allows an unprecedented detailed 
Floquet analysis of the instability, paving the way to rigorous analytical char- 
acterizations of the quasi-collapse unstable manifold. This is a one-dimensional 
manifold, as our numerical results unambiguously confirm. Then, in the regime 
of small time scales, the unstable manifold is fully characterized by a single 
Floquet eigenvector. 




ReX 



Fig. 3. (a) Plot of the Norm, M, of the computed solutions as a function of their 
frequency LVf, for different values of the coupling parameters. The continuations have 
been made starting from the 1-dimensional limit (Ci = 1). For the regions where 
dM/duJb is positive (negative) the continued solutions are stable (unstable). We can 
monitor the change of the linear stability of a solution of a given frequency during 
its continuation in C2 looking at the Floquet spectra. Figures (b) and (c) show the 
Floquet spectra of a discrete breather of frequency Uf, = 3.93 at Ci = 1, C2 = 0.7 
(where dM /dooh < 0) and at Ci = I, C2 = 0.2 (where dM/doJi, > 0), respectively. 



Figure 4 shows the (surface) function Af{uJb, u) for the volume sector of con- 
stant C2{= 0.5). Most noticeably, the critical (threshold) line of bifurcation 
points (^ = 0), as seen in the inset, does not define a monotone function 
ujf^{u). In fact, in the whole interval of < < 1 values, the range of values of 



14 



ul is quite short , indicating the insensitibity of the gross features of the quasi- 
collapse transition to the value of u. However, considering finer details, one 
sees that the threshold curve a;^^(i/) smoothly reaches its slightly larger values 
around midway between the DNLS and the A-L hmits. In other words, inter- 
mediate values of the interpolating (Salerno) parameter u somewhat favour 
the enhancement of the quasi- collapse unstable region. These conclusions are 
in contrast with the stated conclusion (for isotropic lattices) in [22] that the 
Ablowitz-Ladik term increases the stability regime. 




Fig. 4. Surface J\f{u}b, ^) for the case C2 = 0.5. The inset shows the curve v{(jJb) 
corresponding to dM/dub = 0. This curve gives the transition points where the 
discrete breather changes its stabihty character. 

When instability is allowed to develop beyond the fixed point tangent space 
into the nonlinear realm of perturbations, the trajectory obtained by direct 
integration of the equations of motion invariably ends after a transient (of 
time scale given by the real Floquet exponent larger than 1) in a localized 
solution with complex dynamics, the pulson states. These states were char- 
acterized in [3] in the following terms "... where the peak intensity 
oscillates between the central site and its four nearest neighbours... it is not 
known whether these pulson states represent true quasiperiodic solutions to 
the DNLS equation". What makes these trajectories on the unstable nonlin- 
ear quasi-collapsc manifold of much practical relevance and interest is their 
ubiquity: They appear as persistent localized states in the hamiltonian dy- 
namical evolution from a wide variety of initial conditions. Their description 
requires at least two frequencies, namely the internal (genuine breather-like 
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frequency) and the frequency of the oscillations of the breather width around 
a mean width value, which turns out to be less than the width of the unstable 
exact discrete breather. Second and outer shells of neighbours (in both lattice 
axes) also participate in the width oscillations. 
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Fig. 5. (a) Time evolution of the amplitude |<^*n,m| for the localization center 
(n = 11, m = 0) and two adjacent sites (n = 10, m = 0) and (n = 11, m = 1). 
Each amplitude is normalized to its initial value, so that it can be seen how the 
quasi-collapse instability is developed. The parameters of equation (13) are Ci = 1, 
C2 = 0.5, u = 1.0 {fi = 0) and the frequency of the pinned breather is LJb = 3.50. 
When the instability is fully developed, we analyse the final state by means of the 
Power Spectrum S{uj) of the time evolution of 5R[$io,o(i)] (the real part of local- 
ization center). As can be observed in (b) the internal frequency of the breather 
(highest peak in the spectrum) shifts to a higher value (a;^ = 4.03) and the other 
peaks are located at the frequencies of the harmonics resulting from the combina- 
tion of the internal frequency with the frequency {ujqc = 0.78) asociated with the 
amplitude |$n,m| oscillations shown in (a). 



Though a more detailed characterization of the pulson states will be pre- 
sented elsewhere, it is illustrative to consider (figure 5) the power spectrum of 
the field at the central site of a typical trajectory on the unstable nonlinear 
manifold of a quasi-collapsing pinned discrete breather. This shows peaks at 
the combinations ul + jujgc (j = 0,±1,±2...), where cUgc is the frequency of 
the width oscillations characterizing the pulson state, while uj^ > is a fre- 
quency higher than the (initial condition) fixed point frequency Ub- The new 
frequency turns out to be very close to the breather frequency of the same 
(initial) norm on the stable branch. In other words, the instability drives a 
shift of breathing frequency towards the stable branch, while the excess energy 
is transferred to the oscillatory motion of the observable width. This behav- 
ior seems to be the essence of the physical characterization of the nonlinear 
quasi-collapse manifold dynamics. 

The numerical observation of a two-frequency power spectrum for a typical 
pulson state points towards an eventual positive answer to the question (on 
true quasiperiodicity) arised in [3], though, for more details we have to refer to 
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the PhD dissertation [46] . This point serves to illustrate how the high accuracy 
of the fixed point numerical solution provides detailed clues on many still 
unsolved (from a mathematical and physical point of view) questions on two- 
dimensional Schrodinger localization, which are of prospective experimental 
interest in nonlinear (photonic, Josephson, ...) physics technologies. 



4 Breathers moving along the strong coupHng direction. 

Early and current attempts to explore straightaway discrete breather mobility 
in isotropic 2D Schrodinger lattices seem to agree [47] that "kicking" proce- 
dures meet huge difficulties in delivering good mobile solutions, contrary to 
the numerical experiences in ID lattices. We note here that the formal basis 
for those methods [48] takes advantage of the Floquct spectra analysis of exact 
pinned breathers, where the so-called depinning (symmetry-breaking) mode is 
identified. This allows, provided Peierls-Nabarro barriers are small enough, to 
obtain nice numerical ID mobile discrete breathers, by computing trajecto- 
ries from perturbations of the exact pinned breather along the tangent space 
direction specified by the depinning eigenvector. The presence of symmetry- 
breaking instabilities leading to exchange of stability between one-site and 
two-site centered pinned breathers [45,49] and the associated lowering of the 
Peierls-Nabarro barriers to breather displacements, hugely facilitating the suc- 
cess of these procedures when applied to (both hamiltonian and dissipative) 
one-dimensional lattices. [50,51,52] 

In contrast, our "anisotropic lattice" continuation approach takes advantage of 
the availability of exact ID mobile solutions by monitoring the parameter C2 of 
transversal coupling, and then does not rely on how easily one promotes clean 
mobility from pinned localization. In this way we obtain accurate numerical 
{Px,Py = 0, q) fixed points, that is Schrodinger discrete breathers moving along 
the strong coupling direction. In a forthcoming paper we will address the 
(much more difficult) question for arbitrary direction of motion. 

4.1 Structure and stability of (1,0,1) fixed points. 

In figure 6 we visualize the instantaneous real and imaginary components of 
the 2D discrete field profile of a typical (1,0,1) Schrodinger breather. Its struc- 
ture can be seen as the natural extension to two-dimensional lattices of the 
structure of mobile Schrodinger breathers analyzed in [27,28]. The numerical 
solution is spatially asymptotic to a, finely tunec? small-amplitude extended (de- 

localized) radiation state ^m,n (t) when m, n — > 00. The fixed point solution 
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can be thus decomposed as 

^rr,nit) = ^Z'(^) + ^inHt), (31) 



which defines the spatially localized component of the solution. It 

turns out that the spatially delocalized component is a highly localized state 
in the (continuum, in the thermodynamic limit) fc-space of wavevectors. More 

precisely, $m,n (t) is a finite linear combination of (1,0,1) -resonant nonlinear 
{i.e. amplitude-dependent frequency u) 2D planewaves. It can be said that, 
as might be expected, ID Schrodinger breather mobility smoothly persists 
when (strong Ci-coupling) ID chains are coupled transversally. Importantly, 
the numerical continuation for increasing values of the transversal coupling 
C2 proceeds far from the weak coupling regime into where the genuine two- 
dimensional effects start to be manifest, as we will see below. 




Fig. 6. Real, (a) and (b), and imaginary, (c) and (d), of a mobile (1, 0, l)-discrete 
breather of frequency u!/, = 2.712. The parameters of equation (13) are Ci = 1, 
C2 = 0.14 and v = 0.95 (/x = 0.05). The insets in (a) and (c) show the background 
far from the moving core. It can be observed that the wavenumbers in the transversal 
direction are ky = ib7r/2. (b) and (d) show the contour plot for both real and 
imaginary parts. 
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Most noticeable, the SVD-regularized Newton procedure invariably selects the 
values ky = ±7r/2 for all values of C2 and u, and thus the values of for the 
2D resonant planewave are independent of C2 (so it remains equal to the k 
values of the ID (1, 1) fixed point for the uncoupled chain). The appearance of 
an extended background modulation in the transversal direction of ky — ±7r/2 
appears naturally as the best choice to take advantage of approximately ID 
breather propagation along strong coupling direction, for it keeps the value 
of kx favoured by the strong coupling Ci value: Any other value of ky would 
entail a different k^ value. Note however that this provides only a plausibility 
argument for the interpretation of the numerical observation {ky — ±7r/2). 

The high accuracy of the computed solutions allows a detailed analysis of many 
issues concerning 2D Schrodinger breather exact mobility along the strong 
coupling direction. We leave aside in this paper many of them, and focus 
here on how the existence of quasi-coUapse instabilities of pinned Schrodinger 
breathers, for increasing C2-coupling values, influences the stability properties 
of moving (1,0, 1) breathers. In other words, we search here for genuine 2D 
effects on these "strong-coupling -direction" (quasi-lD) moving breathers. 

We have performed an exhaustive exploration of two sectors of the parameter 
space {C2,uJi), u), corresponding to the breather frequency values ujb = 2.5843, 
and Ub = 2.712, by computing the continued (1, 0, 1) fixed point. These values 
of Ub were chosen low enough to allow the analysis of pinned breather quasi- 
coUapse effects on mobility, which occurs at relatively low values of C^^^y) for 
these values of u . 

The Floquet analysis of the computed solutions provides the stability dia- 
grams represented in figures 7. Both show no qualitative differences: There 
are two regions in the (6*2,2^) plane where the (1,0, 1) mobile breather is lin- 
early unstable. The figures are not "schematic": Every point of the plane in a 
fine grid of values of C2 and v has been analyzed, i.e. the Floquet spectrum 
of the computed (1, 0, 1) fixed point is scrutinized, as shown in figures 7.b and 
7.d, where the modulus of the Floquet eigenvalues is shown as a function of 
either v (figure 7.b) or C2 (figure 7.d). 

The first unstable region appears at low values of C2 and intermediate to 
high values of the Salerno parameter z/, i.e. it does not occur close to the 
A-L limit. This unstable region is also bounded above in the direction of C2: 
The variation of the modulus of the unstable Floquet eigenvalue versus the 
transversal coupling parameter C2 shows that the mobile breather becomes 
stable again at larger values of C2, before the second instability at even higher 
coupling takes place. An important observation is that the pinned discrete 
breather of the same frequency is linearly stable at the points in this unstable 
region for (1,0,1) mobile breathers. Thus this instability cannot be ascribed 
to pinned quasi-coUapse effects. 
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Fig. 7. Stability diagram and evolution of the modulus of the Floquet eigenvalues for 
two (1, 0, l)-discrete breather of frequencies oJh = 2.584 (a) and (b), and = 2.712 
(c) and (d). The stability diagram (a) and (c) show two regions where the mobile 
DB becomes unstable. For low values of the coupling C2 there is a subset of values 
of V where the breather suffers from fission (see text and figures (8. a) and (8.b). On 
the other hand for higher values of C2 there is a second region (quasi independent of 
v) where the unstable breather yiels a travelling quasi-collapsing state (see text and 
figures (8.c) and (8.d)). The evolution of the modulus of the Floquet eigenvalues 
along different paths v = 0.50 (b) and C2 = 0.17, 0.16, 0.02 (c) is shown. 



The second transition occurs for values of C2 close to, but slightly higher than, 
the values C*^ of the quasi-collapse of the pinned breather of the same fre- 
quency. We had already seen in the previous section that the quasi-collapse 
transition C^^{v) is only very weakly dependent on the value of i/, and note 
that the same is true for this mobile breather bifurcation. These results suggest 
that this second transition is related to quasi-collapsing phenomena. Signifi- 
cantly, the stability of the (1, 0, 1) mobile breather persists for a small interval 
of coupling values above the pinned breather quasi-collapse. This should be 
regarded as natural, for the mobile breather is a different solution. Note in 
figure 7.b that the modulus of the unstable Floquet eigenvalue, in the inte- 
rior of the unstable region, reaches much higher values than those typical for 
the first type of instability, and decreases for larger values of C2, before the 
breather solution ceases to exist and only plane wave solutions are obtained 
by our numerical method. Note that this behaviour of the unstable Floquet 
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eigenvalue also fits well to the main features of the pinned quasi-coUapse in- 
stability strength, as described by the slope dN'/dui,. From now on we will 
refer to this instability of mobile breathers as the quasi-coUapse instability. 

In the next subsection we characterize both generic types of instability, by 
looking at the details of the unstable manifold associated with each type. As 
we will see, pulson states turn out to play a role in the description of typical 
trajectories on the unstable nonlinear manifolds. 

4.2 Unstable manifold behaviour and ubiquity of pulson states. 

First, we analyze the quasi-collapse instability of (1, 0, 1) mobile breathers. 
The unstable linear subspace in the tangent space of the fixed point is one- 
dimensional. The typical instantaneous profile of the (modulus) unstable Flo- 
quet eigenvector driving the instability is shown in figure 8.d. It is an exponen- 
tially localized 2D profile which decays asymptotically to zero as m,n — > 00, 
i.e. it does not excite radiation. These characteristics are shared by the quasi- 
collapse unstable eigenvector of the pinned breathers, which further reinforce 
the previous considerations leading us to consider this instability as the mobile 
counterpart of the pinned quasi-collapse transition. 

In figure 8.c we have visualized the time evolution of the field modulus con- 
tour plot for a typical trajectory on the unstable manifold. This is obtained by 
direct numerical integration of the equations of motion, from an initial con- 
dition in which a small perturbation along the quasi-collapse eigenvector has 
been added to the unstable fixed point solution. One sees that the breather 
translational motion slows down, and the energy is transferred to width oscilla- 
tions. These oscillations turn out to be more irregular, see figure 9, than those 
observed in section 3 when we inspected typical trajectories on the unstable 
nonlinear manifold of pinned breathers. 

The difference in the character of the width oscillations in both (pinned and 
mobile) cases may be ascribed to the presence of an extended background 
component in the mobile breather solution, which naturally enters into the 
energy transfer taking place during temporal evolution. The slowing down of 
the translational motion continues and eventually the breather pins into a 
convulsive pulson state surrounded by the remaining radiation. 

Now we pay attention to the "low C2" instability of (1, 0, 1) mobile breathers. 
The modulus profile of the unstable Floquet eigenvector that drives this in- 
stability is M-shaped (bimodal), as shown in figure 8.b, and is asymptotic to 
an extended planewave-like profile as m,n 00, i.e. it is not a purely local- 
ized perturbation. It is indeed rather different from the quasi-collapse unstable 
eigenvector analyzed above, which is consistent with the fact that the pinned 
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Fig. 8. Time evolution of two unstable solutions, (a) and (c), of frequency 
omegab = 2.584 and the asociated unstable Floquet eigenvector, (b) and (d) respec- 
tively. Figures (a) and (c) show the time evolution of the contour lines correponding 
to three different values of |$n,m|) in order to visualize the 4-dimensional functions 
|^n,m|(i)- Figures (a) and (b) shows the fission of the breather solution when per- 
turbed along the unstable "M-shaped" Floquet eigenvector plotted in (b). It can be 
seen how a low amplitude pulse emerges and the mobile breather becomes pinned. 
After this transient this low amplitude pulse decays into radiation. The parameter of 
equation (13) are Ci = 1, C2 = 0.08 and v = 0.5 (/i = 0.5). In the case of figures (c) 
and (d) the parameters are the same except for C2 = 0.19. In this case the solution 
is in the "quasi-collapse" unstable region shown in figure (7. a). The final state when 
perturbed along the unstable eigenvector (d) is a travelling breather whose ampli- 
tude oscillates in the same fashion as that of the pinned quasi-collapsing breathers, 
i.e the localization center oscillates out of phase with repect to all the other sites 
on the lattice. 

breather of the same frequency is linearly stable in this region of parameter 
space. As argued above, this instability is not related to quasi-collapse phe- 
nomena, and it does not appear in the region of small values of the Salerno 
parameter z/, close to the A-L limit. 

A typical trajectory on the unstable manifold associated with this instability 
is shown in figure 8. a, where we have plotted the time evolution of the field 
modulus contour plot. We can see there that the mobile breather pins quickly 
while a small pulse moving backwards is ejected, which spreads and finally 
mixes with the remaining delocalized background. However some energy is 
transferred to width oscillations of the pinned breather so that also in this 
case we observe the formation of pulson states surrounded by the remaining 
radiation. As the main difference of this behaviour, with respect to the evo- 
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Fig. 9. Time evolution of the maximum value of the modulus along the cen- 

tral (m = 0) chain and the adjacent (m = 1) one for a mobile (1, 0, 1) breather with 
frequency uJb = 2.584. This magnitude is normalized to the initial value |<l>n,m|(io)- 
Figure(a) shows this evolution for an stable situation (C2 = 0.15, v = 0.5). It can 
be observed how the localization center (m = 0) and is neighbour in the transver- 
sal direction (m = 1) follows two in-phase periodic trajectories in their modulus 
due to the Peierls-Nabarro barrier surpassed during the motion. In contrast, figure 
(b), shows the case when the breather is unstable (C2 = 0.19, v = 0.5). Here the 
quasi-collapse dynamics is manifested while the localization center moves across the 
lattice. As can be observed, the oscillations of the two amplitudes are out of phase 
and the amplitudes of these oscillations are one order of magnitude higher than 
those of figure (a). 



lution observed on the quasi-collapse unstable manifold, is the ejection of the 
small moving pulse, we refer to this instability as fission. 



By increasing the strength of the initial perturbation along the direction of the 
unstable eigenvector, one observes that the size of the ejected pulse increases. 
This observation is consistent with the results reported in [22], where the 
evolution of initial moving gaussian pulses in isotropic 2D Schrodinger lattices 
was studied. These numerical experiences lead the authors to conclude that 
"the characteristic feature of the discrete quasi-collapse of a moving pulse is the 
splitting of the initially moving broad pulse into a track of the standing narrow 
structures ..." (sic). However, we see from our study of the stability of exact 
moving discrete breathers that the fission and the quasi-collapse instabilities 
have different origins and they appear in different regions of parameter space. 
On the other hand, the ubiquitous phenomenon of width oscillations of pinned 
localized structures (pulson states) cannot be ascribed to quasi-collapse. They 
also appear as the preferred way to allocate excess of (localization) energy in 
regions of parameter space far from the quasi-collapse unstable region. 
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5 Conclusions and Prospective Remarks 

We have studied here the dynamics of exact numerical discrete breathers, both 
pinned and mobile, in two-dimensional anisotropic nonlinear Schrodinger lat- 
tices. These solutions are computed using a SVD-regularized Newton method 
by continuation from a set of uncoupled ID chains into increasing non-zero 
values of the coupling in the transversal direction. 

We have performed an extensive exploration in the parameter space (a;^, C2, ly) 
of breather frequency, transversal coupling and Salerno parameter, by comput- 
ing the Floquet spectra of the numerical solutions. We have also computed the 
breather norm J\f{ui„ C2, z^) and further corroborate the general validity of the 
criterion found in [20], namely that the partial derivative dj\f/dub is positive 
for stable pinned breathers. Furthermore, we have analyzed the dynamics on 
the quasi- collapse unstable manifold, where the unstable breather experiences 
a shift in frequency towards the (higher) value of the stable breather with the 
same norm. The excess of energy is coherently transferred to oscillations of 
the breather width, so that the resulting pulson state is characterized by two 
frequencies. 

We have studied discrete breathers moving along the strong coupling direction. 
These solutions are composed of an exponentially localized core on top of an 
extended background which is itself the finite sum of a finite set of nonlinear 2D 
plane waves. The time scales asociated with these plane waves are resonant 
with the core internal frequency as happens in the ID case. In particular, 
the background chooses a finite set of plane waves from a continuous family 
of resonant solutions. The Floquet analysis of these mobile discrete breathers 
reveals the existence of two distinct types of instability. One is the counterpart, 
for mobile breathers, of the quasi-coUapse experienced by pinned breathers. 
The other instability occurs in a region of parameter space where pinned 
breathers are linearly stable. The analysis of the dynamics on the unstable 
manifold show that the excess of energy is partly transferred to a small moving 
pulse, ejected from the center of localization, which justifies the designation 
of a fission instability. However, part of the energy excess is also transferred 
to width oscillations. The appearance of pulson states far from the quasi- 
coUapse regime indicates that the tendency to allocate energy in the form of 
width oscillations is a general 2D feature, not exclusively associated to quasi- 
coUapse instabilities. 

In a future work we will focus on mobility of 2D discrete breathers in an ar- 
bitrary lattice direction. The results obtained here shed light about how this 
mobility can be obtained. In fact, our experiences show that mobility of pinned 
breathers can be induced based on the existence of the extended background 
in the numerically exact mobile solution. On the other hand, the results ob- 
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taincd here and the aforementioned future work may help to design and better 
understand recent numerical experiments reported in [53], concerning the in- 
teraction between high amplitude pinned breathers and mobile ones. These 
experiments provides a possible way for routing and blocking mobile discrete 
breathers via the interaction with the high amplitude pinned ones, resulting 
in a plausible implementation of logical functions. 
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